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1  Introduction 


A  fundamental  understanding  of  the  relationships  between  processing  and  the  overall  quality  and 
ill-service  performance  of  thick-section  thermosetting  composite  structures  is  needed.  The  devel¬ 
opment  of  residual  stresses,  for  example,  is  strongly  influenced  by  processing  history.  Residual 
stresses  can  have  a  significant  effect  on  the  mechanics  and  performance  of  composite  structures  by 
inducing  warpage  or  initiating  matrix  cracks  and  delaminations  [1,2, 3,4, 5]. 

Processing  concerns  associated  with  thermosetting  composites  become  increasingly  important 
for  components  of  appreciable  thickness  [6,7,8].  Perhaps  the  most  familiar  concern  is  an  increase 
in  internal  temperature  resulting  from  the  irrevei  uble  exothermic  chemical  reaction  of  the  ma¬ 
trix  phase.  Liberated  heat  is  slow  to  dissipate  by  conduction  and  may  potentially  raise  internal 
part  temperatures  to  levels  risking  material  degradation.  A  second  concern  relates  to  the  complex 
temperature  and  degree  of  cure  gradients  that  develop  in  thick-sections  during  the  curing  process 
[9,10,11].  These  gradients  induce  non-uniform  curing  within  the  part  that  may  ultimately  lead  to  a 
reduction  in  the  overall  quality  and  in-service  performance  of  the  finished  component.  Non-uniform 
curing  can  result  in  incomplete  consolidation  of  the  part  which,  in  turn,  may  lead  to  ui' desirable 
volume  fraction  gradients  and  entrapped  volatiles  or  voids  [12,13].  Comple.K  cure  gradients  increase 
the  potential  for  process-induced  warpage  and  matrix-microcracking  and  diminished  residual  prop¬ 
erties  of  the  structure. 

The  objective  of  this  study  is  to  gain  a  fundamental  understanding  of  the  curing  process  unique 
to  thick  thermosetting  composites  parts  of  arbitrary  cross-section.  A  two-dimensional  anisotropic 
cure  simulation  analysis  is  presented  which  accounts  for  thermal  and  chemical  interactions  asso¬ 
ciated  with  the  cure.  Several  typical  glass/polyester  and  graphite/epoxy  structural  elements  of 
arbitrary  cross-section  are  analyzed  to  provide  insight  into  the  non-uniform  curing  process  of  thick- 
sections.  Spatial  gradients  in  temperature  and  degree  of  cure,  unique  to  thick-section  composites, 
are  shown  to  be  strongly  dependent  on  part  geometry,  thermal  anisotropy,  the  chemical  cure  kinetics 
and  the  thermal  boundary  conditions  (cure  cycle  and  tooling).  Correlation  with  experimental  mea¬ 
surements  of  through-the-thickness  temperature  profiles  in  glass/polyester  laminates  are  presented 
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for  several  arbitrary  temperature  histories. 


2  Background 

The  autoclave  curing  process  of  thermosetting  composites  has  been  the  subject  of  numerous  in¬ 
vestigations  [6,12,14,15,16,17].  From  early  efforts  to  fabricate  thick-section  composites,  various 
undesirable  effects  were  encountered  that  lead  to  poor  part  quality.  Consequently,  most  studies 
have  sought  to  understand  the  curing  process  on  a  fundamental  level.  Studies  of  the  curing  process 
have  focused  on  the  thermal  and  chemical  interactions,  degree  of  cure  profiles,  viscosity  behav¬ 
ior,  void  formation  and  growth,  and  resin  flow  phenomena  occurring  in  the  composite  under  the 
application  of  a  specified  temperature  and  pressure  cure  cycle  history  [12,13,18].  Other  investiga¬ 
tions  are  more  empirically  oriented,  citing  general  experiences  encountered  in  the  fabrication  and 
manufacture  of  thick-section  thermoset  composite  parts  [7,8].  A  review  of  some  of  the  literature 
associated  with  the  processing  of  thermosetting  composites  is  now  presented.  The  survey  presented 
here  is  not  intended  to  be  exhaustive,  but,  those  authors  cited  are  felt  to  have  made  significant 
contributions  in  the  area. 

Early  investigations  by  Levitsky  and  Shaffer  [9]  focused  on  temperature  and  degree  of  cure 
gradients  that  develop  in  chemically  reacting  isotropic  systems.  Their  one-dimensional  analytical 
solution  with  prescribed  temperature  boundary  conditions  enabled  them  to  investigate  the  influence 
of  various  reaction  kinetic  variables  on  the  curing  process.  They  extended  their  work  to  sho\'  the 
significant  influence  of  temperature  and  degree  of  cure  gradients  on  the  development  of  stress  in 
isotropic  materials  induced  by  non-uniform  curing  [19,20,21]. 

Loos  and  Springer  [12]  developed  a  comprehensive  one-dimensional  simulation  model  to  describe 
tiie  curing  process  of  flat  plate  unidirectional  AS4/3501-6  graphite/epoxy  composite  laminates.  The 
model  integrated  submodels  which  describe  the  fundamental  mechanisms  associated  with  the  cur¬ 
ing  process  such  as  the  thermo-chemical  interactions,  resin  flow  and  void  formation.  Governing 
equations  describing  (.he  curing  process  are  solved  with  an  implicit  finite  difference  method.  Tem¬ 
perature,  degree  of  cure,  resin  flow  and  void  size,  among  other  processing  variables  are  predicted  as 
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a  function  of  the  autoclave  pressure  and  temperature  cure  cycle  history.  Experimental  verification 
of  the  model  was  performed  and  results  are  in  good  agreement  with  simulated  predictions, 

Kays  [6]  lias  conducted  a  comprehensive  three  year  investigation  on  the  processing  issues  unique 
to  large  area  thick-section  laminates.  The  baseline  material  system  was  unidirectional  AS4/3501 
graphite/epoxy.  Cure  simulation  models  were  developed  and  used  in  the  investigation.  Vari¬ 
ous  autoclave  procedures,  cure  monitoring  and  non-destructive  evaluation  (NDE)  techniques  for 
thick-section  laminates  were  developed  and  evaluated.  Contributions  towards  the  development  of 
a  generic  methodology  for  processing  thick-section  laminates  were  made.  Interesting  observations 
reported  in  the  study  were  the  development  of  microcracks  and  delaminations  under  certain  proce.ss- 
ing  conditions,  indicating  the  importance  of  processing  on  the  cure  and  performance  of  thick-section 
composites.  Although  ply-drop  geometries  were  included  in  the  study,  cure  simulation  was  limited 
to  a  one-dimensional  through-the-thickness  analysis. 

Efforts  to  optimize  cure  cycles  for  the  large  scale  manufacture  of  thermosetting  resin  compos¬ 
ites  have  been  attempted,  [14,16,17].  Computer-aided  curing  systems,  utilizing  cure  simulation 
submodels  and  control  feedback  systems,  were  developed.  The  studies  focused  on  reducing  com¬ 
posite  manufacturing  cost  while  improving  part  quality  on  a  reproducible  basis. 

Bogetti  and  Gillespie  [10]  recently  conducted  a  fundamental  study  of  process-induced  residual 
stress  in  thick-section  thermosetting  composites.  A  one-dimensional  cure  simulation  model  is  cou¬ 
pled  to  an  incremental  stress  analysis,  A  constitutive  model  is  proposed  to  describe  '  material 
behavior  during  cure  that  includes  chemical  hardening,  thermal  and  cure  shrinkage  effects.  Residual 
stresses  are  shown  to  be  strongly  influenced  by  gradients  in  temperature  and  degree  of  cure. 

While  cure  simulations  models  for  thermosetting  composites  are  generally  based  on  one-dirnensional 
through-the-thickness  assumptions,  the  cure  of  arbitrarily  shaped  thick-sections  is  significantly  in¬ 
fluenced  by  part  geometry  and  anisotropic  heat  transfer  necessitating  a  two-dimensional  analysis. 
The  present  investigation  contributes  to  the  processing  science  of  thick-section  structural  elements 
encountered  in  the  practical  application  of  thermosetting  composites.  The  two-dimensional  anal¬ 
ysis  developed  is  a  prerequisite  for  studying  the  evolution  of  process-induced  residual  stress  and 
deformation  in  arbitrarily  shaped  thick-sections.  This  work,  therefore,  represents  an  important 
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step  towards  achieving  our  ultimate  goal  of  building  in  quality,  long  life,  predictable  and  reliable 
performance,  durability  and  lower  cycle  costs  of  thick-section  thermosetting  composite  structures 
for  future  Army  systems. 

In  the  following  section  the  problem  formulation  is  discussed.  The  governing  equations  and 
boundary  conditions  describing  the  thermal  and  chemical  interactions  associated  with  the  two- 
dimensional  anisotropic  curing  process  are  presented.  Solution  of  the  pertinent  equations  utilizing 
the  boundary  fitted  coordinate  system  (BFCS)  transformation  technique  in  conjunction  with  the 
alternating  direction  explicit  (ADE)  finite  difference  method  is  then  described. 

3  Analysis 

3.1  Assumptions 

The  analysis  assumes  that  gradients  in  temperature  and  degree  of  cure  normal  to  the  cross-section 
of  the  part  geometry  are  negligible.  This  assumption  is  valid  for  large  area  parts  where  edge  effects 
are  minimal.  All  thermal  and  cure  kinetic  material  pararheters  are  assumed  constant,  independent 
of  time  and  temperature.  It  is  also  assumed  that  no  resin  flow  or  part  thickness  reduction  occurs 
during  the  curing  process.  The  no-resin  flow  assumption  is  reasonable  for  the  large  area,  net  resin 
or  low  bleed  thermosetting  composite  systems  investigated  in  this  study. 

3.2  Heat  Conduction  Equation 

Temperature  solutions  are  based  on  Fourier’s  heat  conduction  equation  for  tv/o-dimensional,  tran¬ 
sient  anisotropic  heat  transfer  with  constant  material  properties  and  an  internal  heat  generation 
source  term.  The  equation  is  well  established  and  is  referenced  here  for  completeness  as  [22]: 


for  T{x,z)  in  D  (1) 

where  D  is  the  domain  of  interest  defined  in  an  orthogonal  (a:,  2)  coordinate  system.  The  term  q 
represents  internal  heat  generation  and  ifc*®,  I'**,  fc®*  are  the  effective  anisotropic  thermal  conduc- 
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tivities,  p  is  the  density,  and  Cp  is  specific  heat  of  the  composite.  T  and  t  are  absolute  temperature 
and  time,  respectively. 

The  coordinate  directions  in  equation  (1)  are  defined  in  a  fixed  (x,  z)  global  coordinate  system. 
Fiber-reinforced  composites  exhibit  anisotropic  thermal  properties  defined  in  a  principle  coordinate 
system  with  coordinate  axes  parallel  and  perpendicular  to  the  fiber  direction.  Fiber  orientation  will 
generally  vary  with  respect  to  the  global  coordinate  system  in  an  arbitrary  shaped  geometry.  The 
effective  anisotropic  thermal  conductivities  in  equation  (1)  are  based  on  the  second  order  tensor 
transformation  of  the  principle  thermal  conductivities  given  by; 

/  kxx  \  f  rni^  mn  \  /  \ 

-mn  *33  (2) 

\  fc®*  /  \  -mn  mn  rri^  —  j  y  *13  j 

In  equation  (2),  m  =:  cos(^),n  =  8in{6)  and  *11,  *33  and  *13  are  the  longitudinal,  transverse 
and  cross-term  thermal  conductivities  of  the  composite  in  its  principle  (1,3)  material  coordinate 
system,  respectively.  Fiber  orientation  within  the  domain  is  assumed  coincident  with  the  curvilinear 
coordinate  system,  (r?,^)i  shown  in  Figure  (1).  In  our  simulations,  the  fiber  orientation  is  coincident 
with  tj  =2  constant  grid  lines  in  the  domain.  This  definition  of  fiber  orientation  requires  that  layers 
witHin  the  laminate  conform  to  the  tool  surface  profile  and  does  not  accommodate  laminate  stacking 
sequence  variations. 

The  rotation  angle,  6,  is  defined  by  the  local  orientation  between  the  fiber  direction  and  the 
global  (x,  z)  coordinate  system.  Transformation  at  a  node  in  the  domain  is  graphically  illustrated 
in  Figure  (1)  where  the  orientation  angle  between  the  principle  (1,3)  and  global  (1,  z)  coordinate 
systems  is  defined. 

3.3  Temperature  Boundary  Conditions 

A  generalized  temperature  boundary  condition  formulation  is  used  to  permit  flexibility  in  the  sim¬ 
ulation  of  the  autoclave  curing  process.  Either  Dirichlet,  Neumann  or  Robin  boundary  conditions 
may  be  enforced  on  the  domain  boundaries.  The  generalized  boundary  condition  is  expressed 


mathematically  as; 


=  0 

/or  T{x,z)  on  dD  (3) 

where  dD  represents  the  domain  boundary  or  part  surface.  The  surface  boundary  temperature  is 
Ts  and  n  is  the  outward  unit  normal  to  the  domain  surface.  The  coefficients  a,  6,  and  c  define 
the  effective  heat  transfer  across  the  domain  boundaries.  The  expression  T{t)  in  equation  (3)  is 
interpreted  as  either  the  ambient  autoclave  temperature  cure  cycle  or  the  actual  boundary  surface 
temperature,  depending  on  the  values  of  a, 6  and  c  specified.  Table  (1)  summarizes  the  three 
possible  boundary  conditions  obtainable  from  this  generalized  formulation. 

The  Neumann  or  insulated  boundary  condition  is  used  when  symmetry  conditions  are  imposed 
to  conserve  computation  time  by  reducing  the  number  of  nodes  required  for  the  simulation.  The 
Dirichlet  or  prescribed  temperature  boundary  condition  is  most  useful  when  experimental  transient 
part  surface  temperatures  are  known.  In  the  presence  of  tooling  or  a  bag  assembly,  these  temper¬ 
ature  profiles  are  not  known  apriori  and  may  be  very  different  from  the  cure  cycle,  In  this  case, 
the  Robin  boundary  condition  is  used  where  {h/k)gjf  defines  the  effective  heat  transfer  coefficient 
quantifying  the  heat  flux  from  the  actual  part  surface  to  the  environment.  This  effective  boundary 
condition  accounts  for  all  thermal  resistance  associated  with  the  tooling  and  bag  assembly. 

3.4  Chemical  Kinetics 

The  term  q  in  equation  (1)  represents  the  instantaneous  heat  generation  per  unit  volume  of  material 
and  is  introduced  to  account  for  the  exothermic  chemical  reaction  associated  with  the  curing  process. 
It  is  a  cure  rate  dependent  term  which  is  evaluated  throughout  the  domain  at  every  time  step  during 
the  incremental  solution. 

The  degree  of  cure,  a,  at  a  material  point  is  defined  as  the  ratio  of  the  cumulative  heat  liberated 
from  the  chemical  reaction,  //(f),  to  the  total  heat  of  the  reaction,  //,.  This  is  mathematically 
expressed  as: 


//(f) 
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(4) 


The  heat  liberated  at  any  point  in  time,  t,  is  expressed  in  integral  form; 

r‘  l,dq. 


(5) 


dq/dt  is  the  rate  of  heat  generation  from  the  cure  reaction.  The  total  b-^at  of  reaction  is  similarly 
expressed  as; 

where  f/  is  the  time  for  complete  reaction.  Heat  generation  at  a  material  point  is  identically  zero 
for  t>  if.  Equations  (4),  (5)  and  (6)  are  combined  to  yield  an  expression  fot  the  cure  rate; 


It  ~  jH^^di’~‘ 


(7) 


The  cure  rate,  da/di,  is  a  time  and  temperature  dependent  function  of  the  reacting  material  system 
and  is  typically  determined  empirically  with  isothermal  Differential  Scanning  Calorimetry  (DSC) 
techniques.  The  cure  rate  is  often  expressed  in  terms  of  time,  temperature  and  degree  of  cure  in 
some  form  of  the  Arrhenius  rate  equation.  Details  of  the  experimental  procedures  for  determining 
cure  rate  expressions  and  the  total  heat  of  reaction  are  found  elsewhere  [23,24].  Once  da/dt  and  H,, 
for  tihe  material  are  known,  the  heat  generation  term  in  equation  (1)  is  straightforwardly  obtained 
by  rearrangement  of  equation  (7); 

fir\ 

(8) 


da 


3.5  Initial  Conditions 


The  model  formulation  is  flexible  enabling  arbitrary  initial  temperature,  Ti,  and  degree  of  cure,  a;, 
distributions  throughout  the  domain  at  the  start  of  the  simulation  to  be  specified.  Uniform  initial 
conditions  are  specified  in  all  simulations  presented  in  this  study; 


T{x,z):zTi  in  D  ai  t  =  0 
Q(a:,2)  =  Qi  m  D  at  i,— 0 

Ti  and  ai  are  taken  to  be  the  ambient  temperature  and  zero,  respectively. 


(9) 
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3.6  Boundary  Fitted  Coordinate  System  Transformation  Technique 

The  boundary  fitted  coordinate  system  (BFCS)  transformation  technique  is  used  to  solve  the  gov¬ 
erning  equation  (1)  and  generalized  boundary  condition  (3).  The  BFCS  technique  is  a  mapping 
technique  in  which  coordinates  in  a  physical  curvilinear  coordinate  system,  (x.z),  arc  transformed 
into  a  computational  rectangular  coordinate  system,  (^,rj).  Motivation  for  using  the  technique 
in  the  present  investigation  is  its  advantage  to  easily  accommodate  complex  shaped  geometries. 
Once  the  governing  equations  and  boundary  conditions  are  transformed  into  the  computational 
domain,  a  straightforward  finite  difference  solution  technique  is  applied.  Solutions  are  obtained  in 
the  computational  domain  and  subsequently  mapped  back  into  the  physical  domain  through  the 
correspondence  of  nodes  in  the  respective  planes.  A  conceptual  representation  of  the  relationship 
between  the  two  coordinate  systems  is  illustrated  in  Figure  (2),  where  a  one-to-one  correspondence 
between  nodes  in  each  coordinate  system  is  noted.  The  transformation  from  the  physical  domain 
into  the  computational  domain  is  based  on  the  Poisson  equation. 

The  BFCS  technique  is  well  documented  [25,26],  so  derivation  and  details  of  the  coordinate 
transformation  procedure  are  omitted  for  brevity.  Utilization  of  the  BFCS  technique  in  this  inves¬ 
tigation  involves  both  a  mesh  generation  procedure  and  coordinate  transformation  of  the  governing 
equation  (1)  and  boundary  condition  (3).  All  the  finite  difference  meshes  presented  in  this  investi¬ 
gation  were  produced  with  the  mesh  generation  code  TGMESH  developed  by  Gilmore  [27]. 


Transformation  of  the  governing  equation  (1)  through  the  Poisson  equation  into  the  computa¬ 
tional  domain  yields: 


d^r  ,  ,  d^T  .  d^T  ,  ,  dT  .  dT  dr 

1  +  =  pc— 


(10) 


where  the  spatial  coordinates  in  the  conjputational  domain  are  ^  and  Finite  dilTerence  formulas 
for  the  spatial  and  time  derivatives  appearing  in  equation  (10)  are  presented  in  the  Appendix.  The 
coefficients,  Ai,  are  functions  of  the  effective  thermal  conductivities,  and  kxz,  and  spatial 

derivatives  in  the  physical  domain  and  are  defined  in  the  Appendix. 

Transformation  of  the  generalized  boundary  condition  (3)  is  also  required.  Details  of  the  trans¬ 
formation  procedure  are  presented  elsewhere  [28].  Consequently,  only  the  final  expressions  are 
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A 

Physical  Curvilinear  Coordinate  System 

I 

Jacobian 

Transformation 


Computational  Rectangular  Coordinate  System 


Figure  2;  Conceptual  Representation  of  the  Physical  and  Computational  Domains 
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presented  here.  Neumann  and  Robin  boundary  conditions  are  functions  of  the  temperature  gra¬ 
dient  on  the  domain  boundary.  The  computational  domain  has  four  separate  faces  comprising  its 
entire  surface  as  indicated  in  Figure  (3).  Each  face  or  side  of  the  domain  has  a  unique  expression 
for  the  temperature  gradient  on  the  boundary  surface,  dTafdh,  given  by: 

^  O’*  ”(-<:)  /“c®  (^1) 

^  on  the  h{-T])  face 

the  coefficients  a,  /?  and  7  and  the  Jacobian  coordinate  transformation,  /,  are  functions  of  the 
spatial  derivatives  in  the  physical  domain,  (*,  2),  and  are  defined  in  the  Appendix.  Substituting 
equations  (11)  into  equation  (3)  yields  the  expressions  for  the  transformed  generalized  temperature 
boundary  conditions  on  each  face  of  the  computational  domain: 

+  bT,  +  cT(t)  =:  0  on  the  n(+0  face 

-f  bTt  +  cT{i)  =  0  on  the  n(-fv)  /“c® 

a^25^  +  6T’,-l-cT(t)  =  0  on  the  n(-0  face  (12) 

+  6T’„  +  cT(t)  s=  0  on  the  n(-q)  face 

Transformations  on  temperature,  (i.e.  T,  and  T{t)),  are  not  required  since  scalars  transform  identi¬ 
cally,  Equations  (10)  and  (12)  constitute  the  transformed  heat  transfer  equation  and  temperature 
boundary  conditioins  governing  the  transient  temperatures  distributions  in  the  composite  during  the 
curing  process.  The  solution  to  the  system  of  re8ultii\g  finite  difference  equations  is  now  presented. 


3.7  Alternating  Direction  Explicit  Finite  Difference  Method 

The  Alternating  Direction  Explicit  (ADE)  finite  difference  method  [29]  is  employed  in  the  solution  of 
equations  (10)  and  (12).  The  ADE  method  is  preferred  over  an  implicit  method  because  substantial 
reduction  in  computation  time  is  realized.  The  ADE  method  does  not  require  the  inversion  of  the 
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tfitnperature  coefficient  matrix  encountered  in  the  implicit  approach.  It  is  also  preferred  over  a 
fully  explicit  method  since  it  has  been  shown  to  be  unconditionally  stable  with  time  step  size,  and 
is  less  sensitive  to  computational  errors  [29], 


The  ADE  finite  difference  method  applies  a  fully  explicit  approach  twice  on  the  computational 
domain  for  each  time  step  increment.  One  sweep  through  the  finite  difference  mesh  is  made  in 
a  systematic  forward  direction  where  a  pseudo-temperature  solution,  u,'j,  at  each  interior  node  is 
defined  explicitly  in  terms  of  adjacent  nodes.  For  the  same  time  step,  another  sweep  is  made  in 
exactly  the  reverse  direction  to  obtain  a  second  psuedo  temperature  solution,  Uij.  The  temperature 
at  node  i,  j,  (Tlj),  is  then  computed  at  each  interior  node  for  the  current  time  step  as  the  arithmetic 
average  of  the  two  pseudo-temperature  solutions: 


Ti  j  =  t 


(13) 


The  computational  mesh  of  nodal  dimensions  m  by  n  is  illustrated  in  Figure  (4).  Central 
differencing  formulas  are  employed  at  all  interior  nodes,  while  forward  and  backward  differencing 
formulas  are  used  on  the  boundaries  nodes.  All  finite  differencing  formulas  are  documented  in  the 
Appendbe. 


Following  the  detailed  procedure  described  by  Barakat  and  Clark  [20],  the  expression  for  uij  at 
the  time  step  <  +  A<  for  an  interior  node  during  a  forward  sweep  is  given  by; 

~  (l/((pCp/At)  -I-  Ai  -f-  A2))i,j  [  (pcp/AOi,jtifj  +  ijjj  -t- 

~  “ij  +  +  (A2)i,j(Ui,j+i  -  u'j  +  -f- 

-  u‘-,,y+i  +  “ilw-i)  + 

(/t4/2).',(lii+,,j  -  uKfj)  +  (14) 


The  Uij’s  at  time  stop  f  -f-  A<  on  the  right  hand  side  of  equation  (14)  all  have  previously  been 
computed  during  that  same  time  step,  (i.e.  those  with  subscripts  (t  -  I,  j),  (ij  -  1),  (i  -f  1,  j  -  1) 
and  (i  —  1,  j  -  1)).  Similarly,  ii,'j  at  time  step  I  -f  A<  for  an  interior  node  during  a  brekward  sweep 
is  given  by; 


=  (l/((pcv/Af)  -b  Ai  +  [  (pcv/A0i„,<j  +  <ijj  + 
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(15) 


+  »i-w)  +  -  «lj  +  + 

-  ■'.'-.j)  +  (/l./2)i, )(»,%■  -  »i^.,)l 


where  the  Vij’s  at  time  step  <  +  At  on  the  right  hand  side  of  equation  (15)  are  those  previously 
computed  during  the  same  time  step,  (i.e.  those  with  subscripts  (/  +  +  l))(i  —  lii  +  1) 

and  (i  -  l,j  -  1)). 

Boundary  nodes  are  treated  in  the  usual  fully  explicit  fashion  since  no  time  derivatives  appear  in 
the  transformed  generalized  temperature  boundary  condition  expressions.  An  explicit  representa¬ 
tion  for  the  boundary  node  temperature,  Tj,  in  terms  of  adjacent  nodes  is  obtained  by  substitution 
of  the  appropriate  one-sided  three-point  differencing  formulas  presented  in  the  Appendix  into  the 
generalized  boundary  condition  equations  (12). 

The  heat  generation  term,  qjj,  in  equations  (14)  and  (15),  is  computed  from  equation  (8), 
rewritten  here  as; 


where  (da/di)jj  is  the  instantaneous  cure  rate  at  node  (i,j)  evaluated  at  time  step  t.  Instantaneous 
cure  rates  are  computed  based  on  instantaneous  temperature  and  degree  of  cure  through  an  em¬ 
pirical  rate  expression.  Temperature  distributions  at  time  t  +  Ai  are  obtained  by  solving  the  finite 
difference  equations  (14)  and  (15)  and  substituting  into  (13)  in  conjunction  with  the  appropriate 
explicit  representations  of  the  boundary  temperatures  defined  in  equation  (12). 


3.8  Degree  of  Cure  Calculation 


The  degree  of  cure  over  each  time  step  increment  during  the  simulation  is  based  on  the  degree  of 
cure  at  the  previous  time  step  and  the  instantaneous  cure  rate  at  the  current  time  step.  The  degree 
of  cure  at  node  (>,j)  for  time  step  <  A<  is  computed  from  the  relationship; 


da\ 


f+A! 


At 


(17) 
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The  value  of  a‘j  is  known  since  it  was  computed  at  the  previous  time  step  increment  and 
is  also  known,  from  the  empirically  determined  cure  rate  expression. 

Temperature  and  degree  of  cure  distributions  throughout  the  domain  are  computed  at  each  time 
step  incremeni  as  a  function  of  the  cure  cycle  temperature  history.  A  computer  code  was  written 
to  perform  the  numerical  calculations.  Results  of  this  investigation  are  discussed  in  the  following 
section. 

4  Results  and  Discussion 

The  cure  simulation  analysis  developed  is  used  to  gain  a  fundamental  understanding  of  the  curing 
process  unique  to  thick-section  composites.  A  summary  of  the  input  is  presented  first.  Verification 
of  the  cure  simulation  analysis  is  demonstrated  by  comparison  with  available  exact  solutions  for 
simulated  temperature  and  degree  of  cure  profiles.  Predicted  temperature  profiles  are  compared  to 
experimentally  measured  values  within  several  2.54  cm  thick  glass/polyester  laminates  exposed  to 
various  temperature  histories.  The  influence  of  the  tooling  on  the  curing  process  is  demonstrated 
through  the  generalized  boundary  condition  formulation.  The  influence  of  laminate  thickness  and 
the  temperature  cure  cycle  ramp  on  non-uniform  through-the-thickness  curing  are  illustrated  to 
demonstrate  the  complex  curing  phenomena  unique  to  thick-section  composites.  Several  typical 
glass/polyester  and  graphite/epoxy  structural  elements  of  arbitrary  cross-section  are  analyzed  to 
show  the  strong  dependence  of  temperature  and  degree  of  cure  gradients  on  geometry  and  material 
thermal  anisotropy. 

4.1  Input  Summary 

Required  input  data  for  the  cure  simulation  analysis  includes;  (1)  solution  details  including  the 
boundary  and  initial  conditions,  (2)  the  finite  difference  mesh  representation  of  the  part  geometry, 
(3)  thermal  properties  and  (4)  a  cure  kinetic  description  of  the  composite  material  system. 
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4.1.1  Solution  Details 


The  desired  cure  cycle  thermal  history,  T(t),  is  segmented  into  descrote  time  step  increments,  At. 
Since  the  ADE  method  used  in  this  investigation  is  unconditionally  stable  with  the  time  step  size, 
little  difficulty  in  obtaining  a  converged  solution  was  encountered.  In  general,  time  step  increments 
between  one  and  two  seconds  proved  to  be  sufficiently  small  to  yield  converged  solutions  for  all  the 
cure  simulation  results  presented. 

The  transient  temperature  boundary  conditions  for  the  analysis  are  based  on  the  specified 
temperature  cure  cycle  history  and  the  effective  heat  transfer  coefficients  a,  6  and  c  defined  in  the 
generalized  boundary  condition  discussed  previously.  Initial  conditions  are  defined  in  equation  (9). 

4.1.2  Mesh  Input 

The  part  geometry  is  discretized  into  a  suitable  finite  difference  mesh  consisting  of  m  by  n  nodes, 
The  physical  (*,«)  coordinates  of  each  node  in  the  mesh  are  used  to  evaluate  the  spatial  derivatives 
appearing  in  transformed  governing  equations.  The  finite  difference  meshes  employed  in  this  inves¬ 
tigation  were  generated  with  a  computer  software  code  TGMESH[27]  that  is  based  on  the  BFCS 
technique.  Three  geometries  typical  of  structural  composite  components  are  used  in  this  investiga¬ 
tion;  a  flat  plate,  a  90*^  right  angle  bend  and  a  ply-drop.  The  finite  difference  mesh  representations 
of  these  geometries  are  shown  in  Figure  (5).  The  flat  plate  and  right  angle  bend  geometry  are 
of  constant  thickness  equal  to  2.54  cm.  The  ply-drop  geometry  consists  of  a  1,27  cm  flat  section 
joined  to  a  2,54  cm  flat  section  by  a  sinusoidal  transition  region.  Fiber  orientation  in  each  geom¬ 
etry  is  assumed  coincident  with  the  (rj  =  constant)  grid  lines.  Solution  accuracy  is  influenced  by 
refinement  of  the  finite  difference  mesh.  Mesh  refinement  was  chosen  to  yield  converged  solutions 
in  all  the  cure  simulations  performed  in  this  investigation. 

4.1.3  Thermal  Properties 

The  thermal  properties  used  in  this  investigation  for  the  glass/polyester  and  graphite/epoxy  com¬ 
posites  are  summarized  in  Table  (2).  The  cross-term  thermal  conductivity,  fcia,  in  the  principle 
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(b)  90®  Right  Angle  Bend  Geometry  Mesh 


(c)  Ply- Drop  Geomeiiy  Mesh 


Figure  5:  Cure  Simulation  Mesh  Geometries 


coordinate  system  of  both  materials  is  identically  zero. 


4.1.4  Cure  Kinetic  Data 

The  complete  description  of  the  cure  kinetics  for  the  composite  includes  the  total  heat  of  reaction 
and  a  description  of  the  rate  of  reaction  as  a  function  of  temperature  and  degree  of  cure.  Reaction 
rate  expressions  for  the  glass/polyester  and  graphite/epoxy  material  systems  are  different  in  form 
due  to  the  inherent  differences  in  the  overall  order  of  the  reaction  kinetics. 

The  glass/polyesfcer  composite  consists  of  CYCOM  4102  polyester  resin,  manufactured  by  the 
American  Cyanamid  Corporation,  and  a  woven  roving  E^-glass  plain  weave  fabric  reinforcement 
containing  approximately  6  yarns  per  inch.  The  reaction  rate  expression  for  the  glass/polyester 
system  is  second-order  overall  [24]: 

^  =  A,exp{-£.EjRT)a^^[\  -  (18) 

R  is  the  universal  gas  constant  and  T  is  absolute  temperature.  The  exponents  mo  and  Mo,  the 
pre-exponential  coefficient,  Ac,  the  activation  energy,  AiJo  and  the  total  heat  of  reaction  are  listed 
in  Table  (3), 

The  graphite/epoxy  composite  contains  Hercules  Corporation’s  3501-6  resin,  reinforced  with 
unidirectional  AS4  graphite  fibers,  The  reaction  rate  expression  for  the  graphite/epoxy  system 
follows  a  markedly  different  form  [12]: 

^  =  (fci +ifc2a)(l -a)(0.47-Q)  Jor  (a  <  0.3) 

^  =  *3(1 -a)  for  (a  >0.3)  (19) 

and  fca  are  defined  by  the  Arrhenius  rate  expressions: 

k\  =  A\exp{—AE\fRT) 

fcj  =  A-iCxp[-i^E2/ RT)  (20) 

k'i  =  A-^expi-AE^/RT) 
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Dirichlet  (prescribed) 

a  =  0 

6=1 

c=  -1 

Neumann  (insulated) 

6  =  0 

c  =  0 

Robin  (convective) 

0=1 

uawsm 

c- 

Table  1:  Generalized  Boundary  Condition  Coefficients 


raoinsyi 

kn/kas 

Glass/Polyester 

Graphite/Epoxy 

EElSuil 

1.26 

9.42x10“* 

2.163x10-“ 

4.457x10“'* 

2 

1,5,10 

Table  2:  Thermal  Properties  for  Glass/Polyester  and  Graphite/Epoxy  Composites 


Glass/Polyester  24 

me 

Oc 

i4e[min.-*] 

AEe[J/mol] 

^r[kJ/kg] 

0.524 

1.476 

3.7x10“ 

1.674x10® 

77.5 

Graphite/Epoxy  [12] 

i4i[min.-*] 

/Ijfinin.-*] 

>l3[min.-*] 

A  El  [J/mol] 
AZv'jij/mol] 
A  £’3(1 /mol] 
77r[kJ/kg] 

2.102x10“ 

-2.014x10“ 

1.960x10® 

8.07x10'* 

7.78x10* 

5.66x10'* 

198.9 

Table  3;  Cure  Kinetic  Parameters  for  Glass/Polyester  and  Graphite/Epoxy  Composites 


20 


The  pi'e-exponential  coefficients  Ai,  and  A3,  the  activation  energies,  A^i,  AE'j  and  Afia,  and 
the  total  heat  of  reaction  for  the  graphite/epoxy  composite  are  summarized  in  Table  (3). 

4.2  Model  Verification 

Cure  simulation  results  were  verified  for  accuracy  by  comparing  predicted  temperature  and  degree 
of  cure  profiles  with  available  exact  solutions.  Temperature  predictions  are  based  on  an  analytical 
one-dimensional  transient  temperature  solution  [22],  Temperature  profiles  in  a  glass/polyester 
laminate  of  thickness  ^=2.64  cm,  initially  at  O^C  and  exposed  to  constant  boundary  temperatures 
of  O^C  at  2=:0(bottom  surface)  and  lO^C  at  2=^(top  surface),  were  compared  at  2=t!/4,  2=^/2 
and  2=:3^/4  locations.  The  cure  simulation  mesh  utilized  is  the  flat  plate  geometry  illustrated  in 
Figure  (6),  employing  insulated  boundary  conditions  on  the  sides  to  isolate  one-dimensional  heat 
transfer  effects.  Excellent  agreement  with  the  exact  solution  is  shown  in  Figure  (6),  providing 
confidence  in  the  accuracy  of  the  cure  simulation  analysis. 

Degree  of  cure  profiles  obtained  from  the  analysis  were  compared  to  independent  results  based 
on  the  cure  rate  constitutive  relations.  Degree  of  cure  profiles  in  a  glass/polyester  flat  laminate 
undfr  uniform,  isothermal  conditions  are  compared  in  Figure  (7).  Results  are  in  excellent  agreement 
which  confirms  the  accuracy  of  this  portion  of  the  analysis. 

4.3  Experimental  Correlation 

Experimental  temperature  distributions  within  2.54  cm  thick  flat  plate  glass/polyester  laminates 
processed  under  various  arbitrary  autoclave  temperature  cure  cycle  histories  were  rorrolatoJ  with 
the  cure  simulation  predictions,  Thermocouples  were  embedded  within  the  laminates  at  various 
locations  through  the  thickness  and  temperature  distributions  were  monitored  during  the  autoclave 
cure  cycle.  The  laminates,  15.24  cm  by  15.24  cm  by  2,54  cm  thick  in  dimension,  were  constructed  by 
stacking  together  approximately  42  plys  of  prepreg.  The  assembly  was  placed  on  a  0.635  cm  thick 
aluminum  caul  plate,  topped  with  a  single  layer  of  bleeder  cloth  and  surrounded  by  an  aluminum 
dam  to  prevent  transverse  resin  flow  and  compaction.  The  entire  assembly  was  placed  within  a  vac¬ 
uum  bag,  sealed  to  the  caul  plate  with  tacky  tape  and  drav/n  to  one  atmosphere  of  vacuum  pressure. 
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igure  7;  Verification  of  the  Degree  of  Cure  Solution  in  a  Glass/Polyester  Laminate 
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The  laminates  were  then  placed  into  the  autoclave  and  subjected  to  various  prescribed  tempera¬ 
ture  cure  cycle  histories.  During  the  cure,  temperature  profiles  at  the  various  locations  within  the 
laminate  were  monitored  with  an  IBM  personal  computer  linked  to  a  KEITHLEY  data  aquisition 
system,  Thermocouple  readings  were  sampled  every  30  seconds  to  ensure  accurate  temperature 
measurements  during  points  in  the  cure  cycle  when  large  temperature  gradients  develop. 

The  thermal  properties  and  cure  kinetic  input  data  pre.sented  previously  for  the  glass/polyester 
system  and  the  flat  plate  mesh  geometry  were  used  in  the  correlation.  A  parametric  study  identified 
the  effective  heat  transfer  coefficient,  (h/fc)e//,  to  be  87[l/m]  on  the  top  surface  and  I25[l/m]  on 
the  bottom  surface  of  the  laminate  assembly  in  contact  with  the  aluminum  tool  surface.  In  addition, 
insulated  boundary  conditions  ware  employed  on  the  sides  to  isolate  through  the  thickness  effects. 
Comparison  of  numerical  and  experimental  transient  temperature  profiles  at  the  center  of  the 
laminates  subjected  to  three  different  cure  cycles  are  presented  in  Figures  (8),  (9)  and  (10).  Good 
agreement  in  this  experimental  study  offers  additional  independent  validation  of  the  cure  simulation 
analysis, 

4.4  Boundary  Condition  Effects 

The  presence  of  the  tooling,  bag  assembly  and  the  internal  autoclave  environment  can  profoundly 
influence  the  heat  transfer  to  the  composite  and  may  ultimately  alter  thn  curing  process  within 
the  part.  The  net  effect  is  that  the  surface  temperature  variation  of  the  part  may  be  significantly 
different  than  the  prescribed  cure  cycle  temperature.  While  convective  heat  transfer  duo  to  heat  flow 
within  the  autoclave  is  not  modeled  explicitly,  flexibility  with  the  generalized  boundary  coiulitiun 
formulation  enables  the  influence  of  the  effective  heat  transfer  to  the  part  on  the  curing  process 
to  be  investigated.  In  the  limit  as  (/»/!:)«//  -+  oo,  the  actual  surface  temp(?rature  of  the  part 
approaches  the  prescribed  cure  cycle  temperature,  T{t),  As  ih/tt)eff  — ♦  0,  an  insulated  boundary 
•condition  is  approached.  The  generalized  boundary  condition  formulation  enables  the  influence  of 
various  tools  and  bagging  assemblies  on  the  curing  process  to  be  investigated. 

The  following  example  demonstrates  the  influence  can  have  on  the  curing  process,  A 

2,54  cm  thick  glass/polyester  laminate  was  subjected  to  cure  cycle  temperature  history  indicated 
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in  Figure  (11)  with  specified  [hlk)t]j  coefficients  on  the  top  and  bottom  surfaces  of  100[l/m]  and 
500[l/m],  respectively.  The  flat  plate  mesh  with  insulated  boundary  conditions  on  the  sides  was 
used  to  isolate  through-the-thickness  effects.  Resulting  temperature  profiles  at  the  surface  and 
center  of  the  laminate  are  shown  in  Figure  (12).  The  higher  {h/k)gjj  coefficient  permits  a  more 
rapid  transfer  of  heat  into  the  laminate,  causing  the  laminate  to  heat  up  and  cure  fYustor.  This  also 
allows  more  heat  to  escape  from  the  laminate  during  the  exotherm,  resulting  in  an  overall  lower 
exotherm. 

Corresponding  degree  of  cure  profiles,  shown  in  Figure  (13),  demonstrate  the  significant  influ¬ 
ence  {hfk)gfj  can  have  on  the  curing  process.  Lower  (h/k)ef/  increases  degree  of  cure  gradients 
during  the  exotherm.  The  effective  heat  transfer  across  the  laminate  surface  plays  an  important 
role  in  the  development  of  residual  stress  and  warpage  during  proce.ssing  [10]. 

4.5  Thickness  Effects 

The  thickness  of  the  laminate  will  strongly  influence  the  curing  process.  The  elfijct  of  increasing 
laminate  thickness  on  the  temperature  and  degree  of  cure  proHles  is  examined  in  glass/polyostor 
laminates  between  1.38  cm  and  6.08  cm  thick.  Flat  plate  mesh  geometries  of  similar  nodal  density  to 
the  flat  plate  mesh  illustrated  in  Figure  (5)  wore  employed.  Prescribed  temperature  boundary  con¬ 
ditions  defined  by  the  glass/polyestcr  cure  cycle  temperalnre  history  illustrated  in  Figure  (11)  were 
imposed  on  the  top  and  bottom  laminate  surfaces.  Insulated  boundary  comlitious  were  enforced 
on  the  sides  to  isolate  through-the-thickness  processing  effects.  Predicted  centerline  temperature 
and  degree  of  cure  profiles  are  illustrated  in  Figures  (14)  and  (15),  respectively. 

Temperature  distributions  at  164  minutes  into  the  cycle,  shown  in  Figure  (16),  domoustrato  that 
while  the  2.64  cm  and  tliiimer  laminates  are  exothcrining,  the  6.08  laminate  is  still  lieatiiig  up  at  thi.s 
point  in  the  cycle.  Similar  gradients  in  the  degree  of  cure  distributions  at  164  luinuteN  are  sliown  in 
Figure  (17).  The  degree  of  cure  gradients  aie  most  severe  at  this  point  in  the  2.54  cm  laminate.  I'lie 
interior  of  the  2.54  cm  laminate  is  essentially  cured.  In  contrast,  the  5.08  cm  laminate  is  ossonlially 
uncured  at  the  interior.  At  a  hater  point  in  the  cure  cycle,  the  tliicker  laminate  will  exotherm  and 
its  distributions  will  reverse  shape  similar  to  the  thinner  laminates,  ultimately  develojiing  more 
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Figure  11:  Typical  Glass/Polyester  and  Graphiie/Epoxy  Temperature  Cure  Cycles 
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Figure  13:  Influence  of  {h/k)gff  on  Degree  of  Cure  Profiles  in  a  Glass/Polyester  Laminate 
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Figure  14;  Influence  of  Thickness  on  Centerline  Temperature  Profiles  in  Glass/Polyester  Laminates 
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Figure  15:  Influence  of  Thickness  on  Centerline  Degree  of  Cure  Profiles  in  Glass/Polyester  Lami¬ 
nates 
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severe  gradients.  These  results  demonstrate  the  complex  temperature  and  degree  of  cure  gradients, 
unique  to  thick-sections,  which  develop  during  the  curing  process.  Bogetti  and  Gillespie  [10]  have 
shown  that  these  gradients  have  a  profound  influence  on  the  evolution  of  process-induced  stress 
and  deformation. 

4.6  Cure  Cycle  Temperature  Ramp  Effects 

The  cure  cycle  temperature  ramp  can  strongly  influence  the  temperature  and  degree  of  cure  gra¬ 
dients  that  develop  during  the  cure.  The  effect  of  the  temperature  cure  cycle  ramp  on  the  degree 
of  cure  gradients  is  examined  in  a  2.54  cm  glass/polyester  laminate  exposed  to  various  cure  cycle 
temperature  ramps.  The  flat  plate  mesh  with  specified  temperature  boundary  conditions  on  the 
top  and  bottom  surfaces  and  insulated  boundary  conditions  on  the  sides  is  used  here.  The  laminate 
was  cured  at  80°C  for  180  minutes  before  being  subjected  to  temperature  cure  cycle  ramps  of  0.0, 
0.25,  0.5  and  1.0  "C/min.,  illustrated  in  Figure  (18).  Predicted  values  of  the  degree  of  cure  at 
the  centerline  minus  the  degree  of  cure  on  the  surface  of  the  laminate,  (oo  -  a,),  are  plotted  in 
Figure  (19)  as  a  function  of  cycle  time  for  the  various  temperature  ramps  investigated. 

Witl*  thit.  particular  laminate  thickness,  2.54  cm,  the  isothermal  cure  under  the  0.0[‘’C/min.] 
ramp  reveals  a  steadily  increasing  positive  value  of  (oe  -  cr,),  typical  of  an  internal  exotherm. 
Above  a  critical  temperature  ramp,  the  surface  temperature  rises  relative  to  the  interior  since  heat 
transfer  by  diffusion  is  low.  Consequently,  the  surface  temperature  initiates  the  cure  reaction.  The 
exothermic  reaction  accelerates  the  cure  and  as  a  result  creates  a  solidification  front  that  sweeps 
the  laminate  from  the  outside  towards  the  interior  of  the  laminate.  Curing  from  the  outside  to  the 
inside  may  potentially  entrap  voids  and  volatile  by-products  of  the  chemical  reaction  and  enhance 
warpage  and  residual  stress  development.  Consequently,  the  temperature  ramp  can  significantly 
influence  the  quality  and  in-service  performance  of  the  part. 

4.7  Anisotropic  Curing  in  Complex  Shaped  Composites 

The  anisotropic  curing  process  in  arbitrary  shaped  geometries  is  now  presented  to  build  another 
level  of  complexity  in  the  cure  simulation  that  enables  realistic  composite  structures  to  bo  modeled. 
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Figure  16;  Teniperatute  Distributions  in  Glass/Polyester  Laminates  at  Exotliorin 
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Figure  18:  Autoclave  Temperature  Cure  Cycle  Ramps 
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Figure  19:  Efl'cct  of  the  Temperature  Ramp  on  Non-uniform  Curing  in  a  2.54  cm  Glass/Polyester 
Laminate 
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The  effect  of  anisotropic  heat  transfer  on  the  curing  process  in  a  2.54  cm  thick  90“  right  angle  bond 
graphite/epoxy  structural  element  is  examined.  The  angle  bend  mesh  depicted  in  Figure  (5)  is 
used.  A  specified  boundary  condition  defined  by  the  grnphite/epoxy  cure  cycle  temperature  history 
illustrated  in  Figure  (11)  was  imposed  on  all  four  sides  of  the  part,  ratios  of  1,  5  and  10, 

typical  of  graphite/epoxy  composites,  were  investigated. 

Temperature  contours  are  shown  in  Figure  (20)  at  146  minutes  into  the  cure  cycle,  when  the 
maximum  exotherm  occurs.  The  isotropic  case,  l:u/fc33  =  1,  exhibits  through-the-thickness  tem¬ 
perature  gradients  resembling  a  one-dimensional  curing  process.  In  contrast,  fcii/fc,33  ratios  of  5  and 
10  demonstrate  that  temperature  gradients  at  the  bend  are  most  severe  at  the  inner  radius  and  are 
accentuated  by  increased  thermal  anisotropy.  This  simulation  clearly  illustrates  the  significance  of 
the  anisotropic  two-dimensional  formulation  used  in  our  cure  simulation  analysis, 

The  curing  process  within  a  glass/polyester  00°  right  angle  bend  part  is  also  examined.  A 
kn/k^a  ratio  of  2  is  assumed  in  the  simulation.  Prescribed  temperature  boundary  conditions 
defined  in  Figure  (21)  are  imposed  on  all  four  sides  of  the  part.  This  cure  cycle  is  similar  to  Llio 
glass/polyoster  cycle  defined  in  Figure  (11)  except  the  intermediau.'  temperature  ramp  to  00°C  is 
omitted. 

Temperature  and  degree  of  cure  contours  at  the  time  of  maximum  exotherm  (1(35  minutes  into 
the  cycle)  are  illustrated  in  Figure  (22).  The  complex  temperature  gradients  indicate  "hot”  regions 
in  each  arm  of  the  part  geometry,  moving  inward  towards  the  bend  as  they  follow  a  sweeping  euro 
front.  Tlio  exotherm  here  is  occurring  on  the  outer  regions  of  the  part,  sweeping  toward  the  interior. 
The  degree  of  cure  contour  at  tills  point  in  the  cycle  indicates  an  outside  to  inside  curing  proces. 
Severe  gradients  are  noted  with  the  surface  fully  cured  and  the  interior  uncured.  This  curing  process 
will  inevitably  entrap  any  voids  or  chemical  by-products,  rL'.'iult  in  incomplete  consolidation  and 
potentially  induce  severe  stress  and  deformation  in  the  pnrt  during  the  curing  process, 

A  final  example  examines  the  cure  of  a  thick-section  glass/polyester  ply-drop  structural  olomont 
illustrated  in  Figure  (5).  Proscribed  temperature  boundary  conditions  defined  in  Figure  (21)  are 
imposed  on  the  top  and  bottom  surfaces  while  insulated  conditions  are  imposed  on  the  sides.  A 
kii/kss  of  2  is  u,sed,  Temperature  and  degree  of  cure  contours  at  the  time  of  maximum  exotherm 
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Figure  21;  Coiidoiised  Glnss/Polye,st;er  Cure  Cycle 
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TEMPERATURE  (C) 


Figaro  22:  Anisotropic  Caring  in  a  2.54  cm  ^'liick  Glass/Polyontor  00"  lliglit.  Anglo  Rond 


(174  minutes;  are  presented  in  Figure  (23).  The  results  identify  a  ’’hot”  region  in  the  transition  area 
with  temperatures  exceeding  the  autoclave  by  45°C.  The  minimum  temperature  occurs  within  the 
interior  of  the  thicker  .section.  "T^he  corresponding  degree  of  cure  contour  indicates  severe  gradient.s 
in  the  thicker  section  where  the  exterior  region  is  fully  cured  and  .ne  in ’crier  is  only  12%  cured. 
In  contrast,  the  thinner  .section  of  the  ply-drop  is  completely  cured  at  thi.s  portion  of  the  ^rocc.sy. 
The  cure  front  is  eTectively  sweeping  from  the  thin  section  to  the  thick  section. 

The  examples  presented  in  this  investigation  clearly  demonstrate  the  complex  cure  behavior 
encountered  in  the  manufacture  of  thick-section  thermoset  structural  components.  Complex  tem¬ 
perature  and  degree  of  cure  gradients,  accentuated  by  part  geometry  and  anisotropic,  heat  transfer, 
demonstrate  the  significant  contribution  made  towards  the  processing  science  of  thicx-sectlon  ther- 
mosets.  The  cure  simulation  analysis  developed  here  provides  the  detailed  prerequiste  information 
needed  for  predicting  process-induced  stress  and  deformation  in  thick-sections  of  arbitrary  cro.ss- 
section.  This  achievement  represents  a  critical  step  in  our  ultimate  goal  of  optimizing  the  quality 
and  in-service  performance  of  thermosetting  structural  components  during  the  manufacturing  pro¬ 
cess. 

ti« 

5  Conclusions 

A  two-dimensional  anisotropic  cure  simulation  analysis  was  developed  to  study  the  complex  curing 
process  within  thick  thermosetting  composites  of  arbitrary  cross-section.  Complex  gradients  in 
temperature  and  degree  of  cure  were  predicted  as  a  function  of  the  autoclave  temperature  history 
and  the  influence  of  the  tool  on  the  curing  process  was  demonstrated.  Several  typical  glass/polyester 
and  graphite/epoxy  structural  elements  of  arbitrary  cross-section  were  analyzed  to  provide  insight 
into  the  non-uniforrn  curing  process  unique  to  thick-sections.  Spatial  gradients  in  degree  of  cure 
are  shown  to  be  strongly  dependent  on  part  geometry,  thermal  anisotropy,  cure  kinetics  and  the 
autoclave  temperature  cure  cycle.  These  spatial  gradients  directly  influence  the  quality  and  in- 
service  performance  of  the  finished  component  by  inducing  warpage  and  residual  stress  during  the 
curing  process.  This  work,  therefore,  represents  an  important  step  toward.s  achieving  our  ultimate 
goal  of  building  in  quality,  long  life,  predictable  and  reliable  performance,  durability  and  lower 
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Figure  23:  Anisotropic  Curing  in  a  Thick  Glass/Polyester  Ply-Drop 
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A  Appendix 

A.l  Finite  Difference  Approximations 

Finite  difference  approximations  to  the  partial  derivatives  used  in  the  analysis  are  presented  below. 
In  the  following  formulas,  the  temperature,  T,  and  the  spatial  coordinates,  x  and  are  repre¬ 
sented  symbolically  as  /.  In  addition,  subscripts  r;  and  ^  denote  partial  differentiation  with  re.spcct 
to  the  computational  coordinates  and  the  subscripts  i  and  j  refer  to  the  nodal  location  on  the 
computational  grid.  Both  central  and  one-sided  differencing  expressions  are  utilized. 

A. 1.1  Contra!  DifForeticUig 

Three-point  central  differencing  formulas  are  used  at  the  interior  nodes  of  the  discretized  domain. 
The  first  order  derivatives  are  given  by; 

/<  =  (/i+x,, •-/<-!, i)/2A^  (21) 

and 

ft)  =  (fij+i  ~  (22) 

Second  order  derivatives  are  given  by; 

-  2/.'.;  +  /,-,,i)/(AO'  (23) 

and 

/r,n-(/h;+i-2/i,i  +  /..;-l)/(Ar?)'  (2d) 

with  the  cross-term  derivative; 

h,  =  (/.>i.j+l  -  fi-ij+i  -  +  /.-i,>-.i)/dA^Arj  (2.5) 

A. 1.2  One-Sided  Differencing 

The  three-point  one-sided  differencing  formulas  used  to  represent  the  temperature  gradientfi  on  the 
boundaries  of  the  dom.ain  in  the  generalized  boundary  condition  formulation  are  givun  below  for 
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each  face  of  the  computational  mesh.  On  the  (<;  =  1)  face,  the  normal  derivative  is  approximated 
by  the  forward  differencing  expression: 


-  fi+2,j)/A^  (26) 

and  on  the  =  m)  face  by  the  backward  differencing  expression: 

ihkj  =  (3/i.i  -  -  /,-2  (27) 

On  the  (r;  =  1)  face,  the  normal  derivative  is  approximated  by  the  forward  differencing  oxpre.ssion: 

ifr,kj  =  +  4/iJ+i  -  /i.;+2)/Ar,  (28) 

and  on  the  (r)  =  n)  face  by  the  backward  differencing  expression: 

ifvkj  =  (3/i.i  -  4/i.i-i  -  (29) 

A{  and  Ar;  are  identically  1  in  the  computational  domain. 

A.  1.3  Explicit  Time  Difforoucing  Formula 

The  explicit  finite  difference  approximation  for  the  time  derivative  of  temperature  is  given  l^y: 

A. 2  Boundary-Fitted  Coordinate  System  Transformation  CoefTicients 

The  coefficients  contained  in  the  transformed  governing  equation  and  generalized  boundary  condi¬ 
tion  are  listed  below.  Derivatives  on  the  physical  domain  coordinates,  x  and  s,  with  respect  to  ( 
and  t]  are  presented  above. 

A. 2.1  Governing  Equation  Coeffleients 

The  Ai  coefficients  appearing  in  the  transformed  governing  equation  are  given  by; 

A\  —  ~  "t"  (31) 


52 


Ai  =  +  k„xl)lJ'^  (32) 

2I3  —  ~  2fcj;»iCr)iC^  +  2k)jf{x^Zri  +  Xf)2^)))/i/  (>^'^) 

/U  =  (ifc^rCi  +  k„C2  +  1k^,Oi)lJ'^  (3'1) 

/it  =  {k^^Ci  +  k,,C^  +  2t„CB)/J^  (35) 

with  the  coafFicients  Ci  given  by: 

Cl  =  ~  4"  ~~ 

C'j  —  ~  "¥  ~  (‘^0 

C3  =  "h  {'^'^) 

Ci  =  ^^^^,,  -  Z,,2ii  +  ~  Ip»l)/ 

Cs  -•  -  X,,X^l  +  {XitXiiJ^  -  x\j^)/j  (40) 

Ce  =  -  ^^)Hh)N  ('ll) 

The  corresponding  Jacobians  appearing  in  the  Ci  expressions  are  given  by: 

J  =  4!jr„  -  «{«,,  (42) 

J(,  =  Xi^Zt,  +  (43) 

Jt)  —  X(,i'iti  "h  x^Zfiti  ~  Xtjt/Z^  —  XtjZ^tj  (44) 


A. 2. 2  Boundary  Conditions  CoofHcionts 

The  coeflicients  used  in  the  normal  direction  temperature  gradients  are  given  by: 


a  = 

('15) 

(i  =  z^z„  +  y,,jy^ 

(■10) 

7  ■•=  x]  4’  y'l 

(47) 
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